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Chapter 1 

Finite-element/boundary-element coupling for 
inflatables: effective contact resolution 

T.M. van Opstal* 


1.1 Introduction 

A fluid-structure interaction technique for the simulation of inflatable structures is 
introduced. The presented finite-element/boundary-element (FEBE) technique cou¬ 
ples an isogeometric finite element discretization of a flexible shell structure with 
an isogeometric boundary element discretization of a Stokes fluid. This technique 
was introduced in ifTOl for planar problems and extended to 3D in ifTTIl . One of the 
marked advantages of this approach is the contact mechanism, lubrication, which is 
an inherent attribute of the flow model. Its role in preventing contact was theoreti¬ 
cally substantiated, but only demonstrated in the planar setting. In the present work, 
we demonstrate the effectiveness of the contact mechanism in the isogeometric and 
three-dimensional setting, discussing various aspects pertaining to the accurate res¬ 
olution of the traction responsible for contact prevention. 

An extensive body of literature already covers the interface of the fields of fluid- 
structure interaction and contact mechanics. If the problem under consideration re¬ 
quires that the gap between two contact surfaces may vanish, an interface tracking 
technique may be applied, sometimes locally Sim, e.g. EEIEI]. If however, a 
problem allows for a finite gap to be maintained between the contacting surfaces. 
Arbitrary Lagrangian-Eulerian a and Space-time m techniques can been used to 
compute problems as challenging as the disreefing of parachute clusters nsiiihiini 
and 1000 spheres falling through a tube Q. In this work a finite (but arbitrarily 
small) gap will likewise be maintained in contact regions. 

The target application for the current approach is inflatable structures. These typ¬ 
ically undergo large deformations with ubiquitous self-contact during the inflation 
process, which often starts from a complex, folded initial configuration. Correct sim¬ 
ulation is contingent to the resolution of every single contact mode throughout the 
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process. To tackle this type of problem, contact treatment based on the lubrication 
effect inherent to viscous flow is highly advantageous: 


1. it avoids explicit contact detection, which is becomes prohibitively expensive and 
non-robust for complex geometries such as folded inflatables. 

2. accuracy is not impaired by artihcial terms added to the model for soft contact 
treatment, such as in e.g. 0. 


Correct simulation of the contact mechanics is contiguous to the accurate resolution 
of the lubrication tractions, which increase sinularly as the contact gap closes. Thus, 
after recapulating the mathematical model in {1.2 accurate approximation and so¬ 
lution through isogeometric analysis and adaptive quadrature schemes are covered 
in 0 The resulting methodology is demonstrated in a numerical test case, a de¬ 
flating balloon, in |1.4| Finally, conclusions are drawn in 11.5 


1.2 Mathematical model of an inflatable structure 


The governing equations of the fluid-structure interaction are elaborated in ifTTlI . 
In this section, this model is recapitulated concisely, and we refer to ifTTl for full 
details. An inflatable structure occupies a boundary segment 17 at time f, which, 
together with a fixed inflow boundary segment fin encloses the interior domain Qt of 
the inflatable. This is schematized in fig. The configuration 9 : (0,T) x Ft 
maps a material point in the reference manifold jc S F to its position on the current 
manifold 9{t,x) G F,, and is sometimes abbreviated as 9t{-) := 9{t,-). Throughout, 
entities related to the reference manifold are underlined. 

The flow on both the interior and exterior domains, resp. Qt and are 


governed by the boundary integral formulation introduced in [1.2.1 The motion of 
the inflatable structure itself is governed by a shell formulation, cf. §1.2.2 Finally, 
the coupling of these two subsystems, as well as a condition enforcing compatibility 


between fluid and structure solutions, is treated in ^ 1.2.3 



Fig. 1.1: Schematic geometry of an inflatable structure for the case d = 3. The ref¬ 
erence and current domains are overlayed, as they coincide over fin. 
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1.2.1 Boundary integral formulation of the fluid 


The flow interior and exterior to the flexible structure is described by Stokes flow, 
supplemented by Dirichlet conditions g on the boundary dQt and suitable radiation 
conditions in the far-fleld. It will be more convenient to pull back the formulation 
to the reference configuration, denoting the boundary condition hy g = godQf and 
the traction jump across <3X2; due to the interior and exterior flows by Tj. The weak 
formulation of the direct boundary integral equation becomes 
given g G H^/^( dQ ), find (Tj, Q G x M such that 

=-ffr(|;,r) V( y/, u) G X K (1.1) 

with bilinear forms 


vf) := ((1 +A)7,0/2 + (1 - 

and boundary integral operators 


oTT J df2 

{Kr^)t:=Jr^l 
4;r Jdn 




with i,k G {0,1,2}, r = \x — y\ the Euclidean distance and 

^/:=((f°-Vr0,)x(f'.Vr0,))-n 


( 1 . 2 ) 


the determinant of the Jacobian of the map 0,. In the above, n and f“, a G {0,1} 
denote the normal and tangent vectors to (3X2 . Furthermore, the composition with 
(0, X 0;) serves to transport the kernel of the integral operator to the reference ma¬ 
nifold. 

In ([T^, the Lagrange multiplier constrains the kernel of the boundary-integral 
operator V). This kernel is related to the undefined average pressure level po, similar 
to the Stokes PDE with pure Dirichlet boundary conditions. Furthermore, 0<A < OO 
denotes the viscosity ratio between the fluids interior and exterior to (3X2, where the 
respective limits A —?► 0 and A —>^ oo correspond to the separate interior and exterior 
problems. 

One of the advantages that the Stokes model gives is an automatic mechanism for 
contact prevention, namely, lubrication. It is shown in |[8l Thms. 33,40] that, as two 
smooth material boundaries advance each other, the fluid opposes this motion with 
a repulsive traction of singular strength of 0{h^^), where h is the distance between 
the boundaries. Thus, the fluid formulation has a built in contact prevention mech- 
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anism which avoids the highly nontrivial explicit treatment of self-contact which 
characterizes folded inflatable structures. 

We conclude this section with a number of remarks: 

1. The boundary integral formulation gives a direct relation between the 
Dirichlet and Neumann data at the boundary, mapping velocities imposed by the 
structure to the tractions imposed on the structure, which is precisely the relation 
typically required in a fluid-structure interaction problem. No meshing and ap¬ 
proximation of the interior domain is thus required, although the solution in the 
interior can be reconstructed a posteriori if desired. 

2. Contrary to a discretization based on a volumetric formulation, it is not more ex¬ 
pensive to treat the flow in the exterior domain concurrently, it simply amounts to 
adjusting the ratio A. Moreover, when the conditions in both fluids are the same, 
A = 1 and we see that the problem actually becomes cheaper to approximate, as 
the dual layer operator Kt need not be assembled. 


1.2.2 Parametrization-free Kirchhoff-Love structure 

We now turn to the equation governing the configuration 0, which is assumed to 
be a member of a Bochner space. A Bochner space L^(0,r;F) contains functions 
/ : [0, r] —>■ Y, such that Jq \\f{t, Oily First and second time derivatives are 

denoted (0 and (0 respectively. We set H/Zq (F) and write the structure 

problem as 

given (Xo,Ai) G x L^{r) and F S L^{0,T;L^{r)), find 9 G {Xo + {0 G 

L^{0,T;X) : 0 G L^(0, T;L^{r)}) such that 

{0r,p}x+W'{0,-p)=F{f,p) a.e.tG(0,T),VpGX, (1.3a) 
0o = Ao, (1.3b) 

00 = At. (1.3c) 

In this problem, F is an external load, Xo and Xi are the initial position and ve¬ 
locity, and W is the internal energy of a Kirchhoff-Love shell, given in terms of a 
configuration 9 as 

W{9) = {eij{9)eki{9) + e^Kij{9)Kk,i9)) ds, (1.4) 


where 


e,y(0) := i(fT,.^.-0„,;,.0m;7), 

Ny(0) ’= ILm&nJ ~ 0m;!Wm;;) 
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are the components of the membrane and bending strains; and 
;= v5ij5ki + 5(1 — v){5ik5ji + dii5jk), 

is the constitutive tensor of the Saint Venant-Kirchhoff material law, which is es¬ 
pecially suited to the anticipated small strains and large rotations HI. The tensor 
rijj := 5ij — UjUj denotes the projection into the tangent space of T and fi j denotes 
the y* component of the gamma gradient of component i of a vector function /. 

The structural model introduces two model parameters, e and v, which represent 
the flexural rigidity and Poisson’s ratio, respectively. 

This form of the Kirchhoff-Love energy has been dubbed pammetrization-free 
in HD as it is formulated on the stress-free state F, without reference to some 
parametrization of this surface. The main advantage is conceptual: eliminating the 
superfluous and arbitrary parametrization from the formulation. Implementation as¬ 
pects for this formulation are treated in HD- 


1.2.3 Transmission conditions and full problem 


Standard Dirichlet-Neumann coupling is employed between the fluid and the struc¬ 
ture. The kinematic condition imposes a velocity on the fluid at the interface and 
reads 


9, on F, 

qn on Ijn, 


(1.5a) 


with q': fin —> K the normal inflow velocity. The dynamic condition imposes conti¬ 
nuity of tractions at the fluid-structure interface 


F(f;.) = -CT(/,Ij,-). (1.5b) 

with Jt the determinant of 9t, according to The model parameter (D can be 
interpreted as a coupling strength and arises in the nondimensionalization HD by 
agglomoration of several parameters (such as the Young’s modulus) from the di¬ 
mensional models. 

In addition to the coupling conditions, an auxiliary condition is to be satisfied, 
connected to the incompressibility of the fluid encapsulated by the structure. The 
configuration provided by the structural equation of motion has to conserve 
the volume, i.e., Qt{9t) = 0 with 

QtW = ^ qdsdt^ , (1.6) 

where the first term is the current volume, the second the initial volume and the last 
the total influx until time t. This auxiliary constraint is imposed by the Lagrange 
multiplier method, and it is proven in HD that this Lagrange multiplier can be seen 
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as the total excess pressure po required to uniformly expand the inflatable structure 
to the correct volume. This constant pressure coincides with the undetermined mode 
constrained from the fluid subproblem 

The above discussion leads to the aggregated FSI model, composed of the weak 
form of the interior and exterior fluid boundary integral equations ( |l.l| i, the La- 
grangian equation of motion for the structure ( |1.3| l augmented with the compatibility 
condition and coupled through the transmission conditions ([rg: 

given q G L^{0,T;L^(Fin)) and (Xo,Zl) xL^{r),find (0,,Tj,po,C) € {Xo + 
X) X X K X K such that 

{Qt,p)x +Vk^(0r;p) +po2!(0f;p) +pCf(0r) + Cl(-4Tj,p)x 

+ an (ij ,w)+b,{Tj,v)+b,{\i,Q- Ff, (ft, = F* (^(f, •)«, r) 

a.e.t e (0,r), V(p,i/,p,t)) eX x xMxM, (1.7a) 

0o=Zo, (1.7b) 

0o=Zi, (1.7c) 

where we extend q and ft by 0 on F and Tin, respectively. 


1.3 Approximation and solution 


In treating the discretization and solution of the mathematical model EL, empha¬ 
sis is put on lubrication and contact mechanics. It was mentioned in (1.2.1 that the 
magnitude of the traction forces grow as 0{hr^) with h the gap size between smooth 
surfaces in contact. Thus, an isogeometric matching-mesh discretization, which pro¬ 
vides the required smoothness and geometrical exactness, is presented in §1.3.1 As 
the gap size h vanishes, local contributions to the boundary integrals become more 
singular on account of the 1 jr and 1 dependence of the kernels of V) and F,. The 


accurate evaluation of integrals in such circumstances is treated in (1.3.2 Finally, 


as the fluid-structure coupling is anticipated to be strong when contact occurs, we 


outline the partitioned iterative procedure for resolving the coupled system in (1.3.3 


1.3.1 Subdivision surface approximation spaces 

Two considerations have lead to the choice of approximation space for this problem. 
Firstly, a C^-continuous basis is desired to have a rotation free discretization of 
the shell structure, and facilitates contact treatment through the lubrication effect. 
The theoretical results depend on a smooth representation of the fluid boundary. In 
that case, there is a time at which the gap width is much smaller than the diameter 
of the contact surface, and the fluid response resembles lubrication. Secondly, an 
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unstructured basis is required, to be able to represent complex geometries, such as 
folded inflatables; and to perform local refinement where it is required to resolve 
fine geometrical features and near-singular tractions in near contact. 

One of the natural candidates is subdivision surfaces with the Catmull-Clark 
subdivision scheme 0. Recently, (truncated) hierarchical refinement has been pre¬ 
sented for this basis in Il22ll20l . One of the freely available CAGD packages in which 
subdivision surface geometries can be designed is Blender™, which has been used 
in this context. 

Catmull-Clark subdivision surfaces can be seen as an extension of cubic b-splines 
to unstructured grids composed of quadrilaterals. Defining the valence of a vertex 
as the number of edges terminating there, we observe that interior vertices of a 
Cartesian mesh are of valency four. An interior vertex of valency other than four 
is called an extraordinary vertex (EV). To extend the definition of the spline basis 
near an EV, a subdivision mask is introduced, in this case the Catmull-Clark mask, 
which preserves -continuity almost everywhere. As such a mask only provides 
an implicit definition of the basis functions at specific points, Stam’s approach M 
is used to evaluate the basis functions and their gradients in numerical quadrature. 
This approach involves I levels of virtual refinement, such that the quadrature point 
lies in a sub-element away from the EV, and can be evaluated as a regular basis 
function. 

As the domains of the fluid and strticture subproblems coincide, it is convenient 
to use the descretization of the structure for the fluid as well. Thus, this approach 
can be seen as a matching mesh approach, such that any geometrical approximation 
is avoided. The advantages are twofold. Eirstly, errors due to geometrical approx¬ 
imation, which are significant in the BEM ifTSll . are now entirely eliminated. One 
can imagine especially in the pervading near-contact modes, that inaccuracies in the 
tractions due to geometrical approximation may easily become disastrous, due to the 
0{h^^) relation to the gap size h. Secondly, we completely eliminate the volumetric 
meshing step, which is otherwise a formidable challenge in the light of the complex 
initial geometries and large deformations characterizing inflatable structures. 


1.3.2 Boundary-element approximation of the fluid 

To evaluate the singular integrals on product domains in the boundary element prob¬ 
lem for the fluid, the regularizing transformations described in ini are used. The 
convergence rates for standard Gauss quadrature rules are known for this set of 
transformations and class of integrands, cf. ca. This suggests an adaptive quadra¬ 
ture scheme in which the quadrature order is increased until some heuristic criterion 
is met. Let be the approximation of the integral over element K with quadrature 
order q. In this work, the criterion — Eir<y| < 10^^ is used. 

It is observed that this adaptive procedure converges and that the amount of eval¬ 
uations requiring an order q decreases rapidly with q, cf. fig. |1.2[ selection of a 
scheme of ^ > 10 is approximately 10"*^ times less likely than selection of the scheme 
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g' = 3. In fig. |1.2| it can also be seen how many levels would be required if an order q 
scheme would be evaluated near an extraordinary vertex. Such an evaluation would 
require the (/ — 1)* power of the subdivision matrix, which can be cached once per 
level. Observe that the number of required levels grows only very slowly and that in 
this case 8 such levels are required. 



Fig. 1.2; The blue histogram shows the number of occurrences of a qua drature order 


q from a sample of 10^ element evaluations in the simulation of ^ 1.4 The red line 
shows the required number of levels / in Stam’s algorithm to evaluate an order q 
Gauss rule. 


1.3.3 Partitioned iterative solution 


Monolithic solution of a discretization of is untractable, as this requires lin¬ 
earization of the coupling terms. In fact, this leads to hypersingular integrals, i.e., in¬ 
tegrands with singularities stronger than those in ( |1.7| i. Despite the strongly coupled 
nature of this problem due to lubrication, a partitioned scheme is therefore preferred. 


Thus, the fluid and structure subproblems, the {yf,v) and (p,fl) terms of (1.7i re¬ 


spectively, are solved separately. The required damping for stability of the coupled 
problem is introduced through the dissipative implicit Euler time-integration scheme 
applied to the structure. 

A subiteration step within our partitioned procedure starts with a linear extra¬ 
polation of the initial data (provided by the solution at the previous time interval), 
which serves as a first approximation of the new coupled solution. Within a fluid- 
structure subiteration, a structural solve is performed first, to ensure compatibility 
of the fluid boundary data. The subiteration is considered converged if (1) the norm 
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of the structure residual is below the tolerance before a Newton solve is performed 
and, in addition, (2) the norm of the fluid update is below that same tolerance. 


1.4 Deflation of a balloon 

The prime objective of the numerical experiment presented here is to assess the 
accuracy of the solution method to resolve the complex self-contact that often arises 
in the simulation of inflatable structures. To this end, as in Qol, an inverted problem 
is considered, namely the deflation of a balloon. 

The initial geometry dQ contains 320 elements in Ijn and 832 elements in F. 
Both these subsets contain 4 extraordinary vertices of valence 3. A small random 
perturbation is applied to F away from (9F so that dQ remains connected. The 
random perturbations causes the structure to wrinkle instead of contract uniformly, 
rendering the structural response closer to reality. The initial configuration is plotted 

in flg.[T3] 



Fig. 1.3: Initial configuration of for the deflation of a balloon. 


The model parameters are u = 0; e = 5.77 lO^'^ corresponding to a flexible shell; 
A = 1 corresponding to identical fluids in the interior and exterior; and ® = 10^^. 
The outflow q is constant in time and has a sine-shaped profile with a magnitude 
such that the volume would vanish at T = 2*^. The time step size is T = 4 and the 
tolerance in the partitioned solver TOL is set to 10^®, settings for which only 1-2 
subiterations are required at each time level. 
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(a) t = 4.0 10^, m{Qt)/m{Qo) = 0.80 


(b) t = 8.0 10^ m(n,)/m{Qo) = 0.61 




-3 


1,2 



-3 


1\2 



(c)f = 1.20 10^m(^2,)/m(^2o) =0.41 


(d)f=1.31 10^m(^2,)/m(^2o) =0.26 


Fig. 1.4: Snapshots of the configuration and fluid traction at different time levels 
during the deflation process. The color coding corresponds to the vertical component 
of the traction, Tj 2 and m{A) denotes the volume of set A. 
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In fig. |1.4| snapshots of the deflation process are shown. Recall that the plotted 
traction component Tj 2 does not include the contribution of the total excess pressure 
PQn. The action of the lubrication effect can clearly be observed along horizontal 
folds, where the vertical component is positive (red) on upper surfaces and negative 
(blue) on lower surfaces and close to zero elsewhere in figs. 1.4a- 1.4c Thus, the 


collapse of the structure under the action of the total excess pressure pan, which is 
uniform over 17 , is observed to be counteracted locally by Tj to prevent contact. 

As the limits of the color bars are identical, it is clear that the magnitude of these 
lubrication forces increases as the deflation proceeds. At the same time, the folds 
in the fabric be come gradually sharper. The deflation process is terminated at f = 
1.31 10^, cf. fig. 1.4d As can be seen from this snapshot, a large downward traction 
is exerted on the part of the structure near Tin, where the outflow condition ( |1.5a| l is 
prescribed. The outflux q does not correspond to the velocity of lin. The segments 
therefore seem to recede and no lubrication force is instigated. Ultimately, this self¬ 
intersection is therefore not due to insufficient resolution of the fluid response, but 
rather to the artificial problem setup where a structure is deflated instead of inflated. 


1.5 Conclusions 

In this chapter we have presented a FEBE method for the simulation of inflatable 
structures. A Kirchhoff-Love shell with low flexural rigidity, discretized with the 
finite element method (using a Catmull-Clark subdivision surface basis) is coupled 
with a boundary element method discretization of both the interior and exterior fluid. 
It is advantageous to let the fluid inherit the structural mesh, as this leads to a ge¬ 
ometrically exact matched discretization at the fluid-structure interface, which in 
turn enables accurate approximation of contact forces. The coupled system may be 
advanced in a partitioned iterative fashion. 

We have furthermore assessed the capabilities of this method in the presence of 
multiple modes of self-contact, by considering a deflation problem, in which the 
fluid response is anticipated to be dominated by lubrication forces. It was observed 
that the computed traction forces effectively prevent self-intersection of the struc¬ 
ture, although very large deformations are sustained. 


References 

1. M. Bischoff, W. A. Wall, K.-U. Bletzinger, and E. Ramm. Encyclopedia of Computational 
Mechanics, volume 2 Structures, ehapter 3 Models and Finite Elements for Thin-walled Struc¬ 
tures, pages 59-137. Wiley, 2004. 

2. E. Catmull and J. Clark. Recursively generated b-spline surfaces on arbitrary topological 
meshes. Computer-Aided Design, 10(6):350-355, 1978. 

3. M.-C. Hsu, D. Kamensky, Y. Bazilevs, M.S. Sacks, and T.J.R. Hughes. Fluid-structure inter¬ 
action analysis of bioprosthetic heart valves: significance of arterial wall deformation. Com- 








12 


T.M. van Opstal 


putational Mechanics, 54:1055-1071, 2014. 

4. T.J.R. Hughes, W.K. Liu, and T.K. Zimmermann. Lagrangian-Eulerian finite element for¬ 
mulation for incompressible viscous flows. Computer Methods in Applied Mechanics and 
Engineering, 29:329-349, 1981. 

5. A.A. Johnson and T.E. Tezduyar. Mesh update strategies in parallel finite element computa¬ 
tions of flow problems with moving boundaries and interfaces. Computer Methods in Applied 
Mechanics and Engineering, 119:73-94, 1994. 

6. A.A. Johnson and T.E. Tezduyar. Advanced mesh generation and update methods for 3d flow 
simulations. Comput. Mech., 23:130-143, 1999. 

7. D. Kamensky, M.-C. Hsu, D. Schillinger, J.A. Evans, A. Aggarwal, Y. Bazilevs, M.S. Sacks, 
and T.J.R. Hughes. An immersogeometric variational framework for fluid-structure interac¬ 
tion: Application to bioprosthetic heart valves. Computer Methods in Applied Mechanics and 
Engineering, 284:1005-1053, 2015. 

8. T.M. van Opstal. Numerical methods for inflatables with multiscale geometries. PhD thesis, 
Eindhoven University of Technology, 2013. 

9. T.M. van Opstal and E.H. van Brummelen. A finite-element/boundary-element method for 
large-displacement fluid-structure interaction with potential flow. Computer Methods in Ap¬ 
plied Mechanics and Engineering, 266:57-69, http://dx.doi.Org/10.1016/j.cma.2013.07.009 
2013. 

10. T.M. van Opstal, E.H. van Brummelen, R. de Borst, and M.R. Lewis. A finite- 
element/boundary-element method for large-displacement fluid-structure interaction. Com¬ 
putational Mechanics, 50(6):779-788, 2012. 

11. T.M. van Opstal, E.H. van Brummelen, and G.J. van Zwieten. A flnite-element/boundary- 
element method for three-dimensional, large-displacement fluid-structure interaction. Com¬ 
puter Methods in Applied Mechanics and Engineering, 284:637-663, 2015. 

12. S.A. Sauter and C. Schwab. Boundary Element Methods. Springer Series in Computational 
Mechanics. Springer-Verlag, 2011. 

13. M.A. Scott, R.N. Simpson, J.A. Evans, S. Lipton, S.P.A. Bordas, T.J.R. Hughes, and T.W. 
Sederberg. Isogeometric boundary element analysis using unstructured T-splines. Computer 
Methods in Applied Mechanics and Engineering, 254:197-221, 2013. 

14. J. Stam. Exact evaluation of Catmull-Clark subdivision surfaces at arbitrary parameter values. 
In proceedings of the 25‘^ annual conference on computer graphics and interactive techniques, 
pages 295^04. ACM, 1998. 

15. K. Takizawa, M. Fritze, D. Montes, T. Spielman, and T.E. Tezduyar. Fluid-structure interaction 
modeling of ringsail parachutes with disreefing and modified geometric porosity. Computa¬ 
tional Mechanics, 50:835-854, 2012. 

16. K. Takizawa, T. Spielman, and T.E. Tezduyar. Space-time FSl modeling and dynamical anal¬ 
ysis of spacecraft parachutes and parachute clusters. Computational Mechanics, 48:345-364, 
2011 . 

17. K. Takizawa, T.E. Tezduyar, C. Boswell, R. Kolesar, and K. Montel. FSl modeling of the 
reefed stages and disreefing of the Orion spacecraft parachutes. Computational Mechanics, 
54:1203-1220, 2014. 

18. T.E. Tezduyar. Stabilized finite element formulations for incompressible flow computations. 
Advances in Applied Mechanics, 38:1-44, 1992. 

19. T.E. Tezduyar, K. Takizawa, C. Moorman, S. Wright, and J. Christopher. Space-time finite 
element computation of complex fluid-structure interactions. International Journal for Nu¬ 
merical Methods in Fluids, 64:1201-1218, 2010. 

20. X. Wei, Y. Zhang, T.J.R. Hughes, and M.A. Scott. Truncated hierarchical catmull-clark sub¬ 
division with local refinement. Computer Methods in Applied Mechanics and Engineering, 
291:1-20, 2015. 

21. T. Wick. Flapping and contact FSl computations with the fluid-solid interface- 
tracking/interface-capturing technique and mesh adaptivity. Computational Mechanics, page 
published online, 2013. 

22. U. Zore, B. Jlittler, and J. Kosinka. On the linear independence of (truncated) hierarchical 
subdivision splines. Geometry-v Simulation Report, 17, 2014. 



